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Progress in the Effective Field Theory of two and three nucleon systems is sketched, 
concentrating on the low energy version in which pions are integrated out as exphcit 
degrees of freedom. Examples given are calculations of deuteron Compton scattering, 
lO ■ three body forces and the triton, and nd partial waves. 

O . 1. Foundations of Effective Nuclear Theory 

^""^ This presentation is a cartoon of the Effective Field Theory (EFT) of two and three 

nucleon systems as it emerged in the last three years, using a lot of words and figures, and 
a few cheats. For details, I refer to the bibliographies in [|I| as well as a recent review , 
^ I and papers with J.-W. Chen, R.P. Springer and M.J. Savage 0, P.F. Bedaque 0,^, and 
fn ■ F. Gabbiani |^. I will mainly concentrate on the theory in which pions are integrated 
out as explicit degrees of freedom, but also comment on the extensions to include pions. 
K> ! B. Holstein's talk provides an introduction to EFT, and E. Epelbaum's contribution 
H ' investigates one of the two proposals to include pions in more detail. 

Effective Field Theory methods are largely used in many branches of physics where a 
separation of scales exists 0. In low energy nuclear systems, the scales are, on one side, 
the low scales of the typical momentum of the process considered and the pion mass, 
and on the other side the higher scales associated with chiral symmetry and confinement. 
This separation of scales produces a low energy expansion, resulting in a description of 
strongly interacting particles which is systematic and rigorous. It is also model indepen- 
dent (meaning, independent of assumptions about the non-perturbative QCD dynamics): 
Given that QCD is the theory of strong interactions and that chiral symmetry is broken 
via the Goldstone mechanism, Wilson's renormalisation group arguments show that there 
is only one local low energy field theory which originates from it: Chiral Perturbation 
Theory and its extension to the many-nucleon system discussed here. 

1.1. The Lagrangean 

Three main ingredients enter the construction of an EFT: The Lagrangean, the power 
counting and a regularisation scheme. First, the relevant degrees of freedom have to be 
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identified. In his original suggestion how to extend EFT methods to systems containing 
two or more nucleons, Weinberg [Q noticed that below the A production scale, only 
nucleons and pions need to be retained as the infrared relevant degrees of freedom of low 
energy QCD. Because at these scales the momenta of the nucleons are small compared 
to their rest mass, the theory becomes non-relativistic at leading order in the velocity 
expansion, with relativistic corrections systematically included at higher orders. The 
most general chirally invariant Lagrangean consists hence of contact interactions between 
non-relativistic nucleons, and between nucleons and pions, with the first terms of the form 

£nn = N\ido + —)N+ ^tr[(9^St)(9'^S)] + gAN^A-aN - (1) 
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where A^ = (^) is the nucleon doublet of two-component spinors and P* is the projector 
onto the iso-scalar-vector channel, Paa^ = 7f(c"2C"*)a(''"2)a- The iso-vector-scalar part 
of the NN Lagrangean introduces more constants Ci and interactions and has not been 
displayed for convenience. The field ^(x) = y/T, = e}^l^^ describes the pion. D^ = 9^ + V^ 
is the chirally covariant derivative, and A^ = |(^<9^^^— ^^<9^0 (^) '^^ axial (vector) pionic 
current. The interactions involving pions are severely restricted by chiral invariance. As 
such, the theory is an extension of Chiral Perturbation Theory and Heavy Baryon Chiral 
Perturbation Theory to the many nucleon system. Like in its cousins, the coefficients of 
the low energy Lagrangean encode all short distance physics - branes and strings, quarks 
and gluons, resonances like the A or p - as strengths of the point-like interactions between 
particles. As it is not possible yet to derive these constants by solving QCD e.g. on the 
lattice, the most practical way to determine them is by fitting to experiment. 

1.2. The Power Counting 

Because the Lagrangean (|^) consists of infinitely many terms only restricted by symme- 
try, an EFT may at first sight suffer from lack of predictive power. Indeed, as the second 
part of an EFT formulation, predictive power is ensured by establishing a power count- 
ing scheme, i.e. a way to determine at which order in a momentum expansion different 
contributions will appear, and keeping only and all the terms up to a given order. The 
dimensionless, small parameter on which the expansion is based is the typical momentum 
Q of the process in units of the scale Ann at which the theory is expected to break down. 
Values for Ann and Q have to be determined from comparison to experiments and are 
a priori unknown. Assuming that all contributions are of natural size, i.e. ordered by 
powers of Q, the systematic power counting ensures that the sum of all terms left out 
when calculating to a certain order in Q is smaller than the last order retained, allowing 
for an error estimate of the final result. 

For extremely small momenta, the pion does not enter as explicit degree of freedom, 
and all its effects are absorbed into the coefficients Ci. This Effective Nuclear Theory 
with pions integrated out (ENT(^)) was recently pushed to very high orders in the two- 
nucleon sector P] where accuracies of the order of 1% were obtained. It can be viewed 
as a systematisation of Effective Range Theory with the inclusion of relativistic and 
short distance effects traditionally left out in that approach. Because of non-analytic 



contributions from the pion cut, the breakdown scale of this theory must be of the order 
'^Ann ~ 1^77- For the ENT with exphcit pions, we would suspect the breakdown scale to 
be of the order of M^ — M, since the A is not an explicit degree of freedom in (|l|). 

Even if calculations of nuclear properties were possible starting from the underlying 
QCD Lagrangean, EFT simplifies the problem considerably by factorising it into a long 
distance part which contains the infrared-relevant physics and is dealt with by EFT 
methods and a short distance part, subsumed into the coefficients of the Lagrangean. 
QCD therefore "only" has to provide these constants, avoiding full-scale calculations of 
e.g. bound state properties of two nuclear systems using quarks and gluons. EFT provides 
an answer of finite accuracy because higher order corrections are systematically calculable 
and suppressed in powers of Q. Hence, the power counting allows for an error estimate of 
the final result, with the natural size of all neglected terms known to be of higher order. 
Relativistic effects, chiral dynamics and external currents are included systematically, and 
extensions to include e.g. parity violating effects are straightforward. Gauged interactions 
and exchange currents are unambiguous. Results obtained with EFT are easily dissected 
for the relative importance of the various terms. Because only S'-matrix elements between 
on-shell states are observables, ambiguities nesting in "off-shell effects" are absent^. On 
the other hand, because only symmetry considerations enter the construction of the La- 
grangean, EFTs are less restrictive as no assumption about the underlying QCD dynamics 
is incorporated. Hence the proverbial quib that "EFT parameterises our ignorance" . 

In systems involving two or more nucleons, establishing a power counting is complicated 
because unnaturally large scales have to be accommodated: Given that the typical low 
energy scale in the problem should be the mass of the pion as the lightest particle emerging 
from QCD, fine tuning is required to produce the large scattering lengths in the S wave 
channels (1/a ^° = —8.3 MeV, 1/a ^^ = 36 MeV). Since there is a bound state in the ^Si 
channel with a binding energy B = 2.225 MeV and hence a typical binding momentum 
7 = yMB ~ 46 MeV well below the scale Ann at which the theory should break down, 
it is also clear that at least some processes have to be treated non-perturbatively in order 
to accommodate the deuteron, i.e. an infinite number of diagrams has to be summed. 

For simplicity, let us first turn to the EFT in which the pion is integrated out into the 
coefficients '^Ci of the Lagrangean. Here, a way to incorporate this fine tuning into the 
power counting was suggested by Kaplan, Savage and Wise [|12|, and by van Kolck [jl3|. 



At very low momenta, contact interactions with several derivatives - like p^ '^Ci - should 
become unimportant, and we are left only with the contact interactions proportional to 
^Cq. The leading order contribution to nucleons scattering in an S wave comes hence 
from four nucleon contact interactions and is summed geometrically as in Fig. |l| in order 
to produce the shallow real bound state, i.e. the deuteron. 



^Here, an extended note is in order. It has been noticed twenty years ago that an experiment cannot 
measure ofF-shell effects because in a field theory, no ambiguities arise from how we choose to define 
the interpolating fields lin]. We can consistently neglect all operators that vanish by the equations of 
motion because each of them makes a contribution to an observable that has exactly the same form as 
higher dimension operators that are present in the theory. Further, such operators can be removed by 
field redefinitions, and therefore it is consistent to work with a Lagrange density that does not contain 
operators that vanish by the equations of motion. Explicit examples can be found e.g. in [^ and (ll| 
App. B]. Of course, particles in loops still propagate off-shell in an EFT, just as usual. 





-'Co 



i-<o 

Figure 1. Re-summation of the contact interactions into the deuteron propagator. 

How to justify this? Dimensional analysis allows any diagram to be estimated by scaling 
momenta by a factor of Q and non-relativistic kinetic energies by a factor of Q'^/M. The 
remaining integral includes no dimensions and is taken to be of the order Q^ and of 
natural size. This scaling implies the rule that nucleon propagators contribute one power 
of M/Q"^ and each loop a power of Q^/M. Assuming that 

fCn ~ ^Co (2) 

the diagrams contributing at leading order to the deuteron propagator are indeed an infi- 
nite number as shown in Fig. |I], each one of the order 1/{MQ). The deuteron propagator 

(3) 



has the correct pole position and cut structure when one chooses ^Co{fi) = ^ -^. 

^Cq becomes dependent on an arbitrary scale /i because of the regulator dependent, 
linear UV divergence in each of the bubble diagrams. Indeed, when choosing /i ~ Q, the 
leading order contact interaction scales as demanded in (^. As expected for a physical 
observable, the NN scattering amplitude becomes independent of //, the renormalisation 
scale or cut-off chosen. All other coefficients ^Ci can be shown to be higher order, so that 
the scheme is self-consistent. Observables are independent of the cut-off chosen. 

The linear divergence does not show in dimensional regularisation as a pole in 4 dimen- 
sions, but it does appear as a pole in 3 dimensions which we subtract following the Power 
Divergence Subtraction scheme [l^. Dimensional regularisation is chosen to explicitly 



preserve the systematic power counting as well as all symmetries (esp. chiral invariance) 
at each order in every step of the calculation. At leading (LO), next-to-leading order 
(NLO) and often even N3L0 in the two nucleon system, it also allows for simple, closed 
answers whose analytic structure is readily asserted. Power Divergence Subtraction moves 
hence a somewhat arbitrary amount of the short distance contributions from loops to 
counterterms and makes precise cancellations manifest which arise from fine tuning IK 



1.3. Including Pions 

The power counting of the zero and one nucleon sector of the theory are fixed by Chiral 
Perturbation Theory and its extension to the one baryon sector. Therefore, the question 
to be posed is: How does the power counting of the contact terms ^C^ change above the 
pion cut, i.e. when the pion must be included as explicit degree of freedom? We do not 
know yet, but the following two promising proposals are on the market. 

If pulling out pion effects does not affect the running of Co too much, one surprising 
result arises: Because chiral symmetry implies a derivative coupling of the pion to the 
nucleon at leading order, the instantaneous one pion exchange scales as Q^ and is smaller 



than the contact piece ^^^Cq ~ ^Cq ~ Q^^. Pion exchange and higher derivative contact 
terms appear hence only as perturbations at higher orders. The LO contribution in this 
scheme is still given by the geometric series in Fig. |l|. In contradistinction to iterative 
potential model approaches, each higher order contribution is inserted only once. In this 
scheme, the only non-perturbative physics responsible for nuclear binding is extremely 
simple, and the more complicated pion contributions are at each order given by a finite 
number of diagrams. For example, the NLO contributions to the deuteron are the one 
instantaneous pion exchange and the four nucleon interaction with two derivatives. Fig. ^. 
The constants are determined e.g. by demanding the correct deuteron pole position and 
residue. This approach is known as the "KSW" counting of ENT, ENT(KSW) llil 




Figure 2. The next-to- leading order in ENT (KSW). 



If on the other hand, the power counting for '^Cq is dramatically modified when one 
includes pions, both the Cq interactions and the one pion exchange might have to be 
iterated. In his paper, Weinberg |Q therefore suggested to power count not the amplitude 
- as is usually done in EFT - but the potential, and then to solve a Schrodinger equation 
with a chiral potential, as pictorially represented in Fig. |^. E. Epelbaum's talk gives 
more details on the results obtained so far in this approach. 
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Figure 3. The leading order in ENT (Weinberg). 

Both power countings are presently under investigation for consistency and convergence, 
and each has its advantages and shortcomings 0]. If both approaches are self-consistent. 
Nature will tell us which it chose. Although in general process dependent, the expansion 
parameter is found to be of the order of | in most applications, so that NLO calculations 
can be expected to be accurate to about 10%, and N2L0 calculations to about 4%. In 
all cases, experimental agreement is within the estimated theoretical uncertainties, and 
in some cases, previously unknown counterterms could be determined. 



2. Applications 

2.1. Deuteron Compton Scattering 

The first example we turn to demonstrates not only how simple it is to compute pro- 
cesses involving both external gauge and exchange currents in a self-consistent way, but 
especially how powerful the power counting is in estimates of theoretical uncertainties. 



Using the power counting ENT(KSW) in which pions are perturbative, the elastic 
deuteron Compton scattering cross section [Q to NLO is parameter-free with an accu- 
racy of 10%. Contributions at NLO include the pion graphs that dominate the electric 
polarisability of the nucleon aE from their — behaviour in the chiral limit. Because 
they are NLO, the power counting predicts their contribution to be sizeable, namely 20 
to 30%. The comparison with experiment in Fig. ^ shows good agreement and therefore 
confirms the HBxPT value for aE- The deuteron scalar and tensor electric and magnetic 
polarisabilities are also easily extracted §]. 
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Figure 4. The differential cross section for elastic 7-deuteron Compton scattering at inci- 
dent photon energies of E^ = 49 MeV and 69 MeV in EFT(KSW) [^, no free parameters. 
Dashed: LO; long dashed: NLO without the graphs that contribute to the nucleon polaris- 
ability; solid curve: complete NLO result. The estimate of the accuracy of the calculation 
at NLO (±10%) is indicated by the shaded area. 



Compton scattering has also been investigated in ENT (Weinberg) [^ and in ENT(^) 
[0. Surprisingly, the results all agree within the theoretical uncertainty of the calcu- 
lations even at relatively high momenta. This seems to suggest that non-analytic pion 
contributions from meson exchange diagrams are small. 

2.2. Three Body Forces and the Triton 

Since all interactions permitted by the symmetries must be included into the La- 
grangean, EFT dictates that in the three body sector, interactions like 



Abody = H {N^Nf 



(4) 



with unknown strength H are present. Turning again to ENT(^), we therefore have to 
ask at which order in the power counting they will start to contribute. Bedaque, Hammer 
and van Kolck |]TB[ found the surprising result that an unusual renormalisation makes the 
three body force of leading order in the triton channel. 

In order to understand this finding, let us first consider the LO diagrams which come 
from two nucleon interactions. The absence of Coulomb interactions in the nd system 
ensures that only properties of the strong interactions are probed. All graphs involving 



only '^Cq interactions are of the same order and form a double series which cannot be 
written down in closed form. Summing all "bubble-chain" sub-graphs into the deuteron 
propagator, one can however obtain the solution numerically from the integral equation 
pictorially shown in Fig. |^ within seconds on a personal computer. 
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Figure 5. The double infinite series of LO "pinball" diagrams, some of which are shown 
in the first line, is equivalent to the solution of the Faddeev equation of the second line. 

Because one nucleon is exchanged in the intermediate state, each diagram is of the 
order of the nucleon propagator, i.e. Q~^, while the first three buddy force (|) seems to 
be of order Q". We are therefore tempted to assume that three body forces are at worst 
N2L0, i.e. of the order 10%. In the quartet channel, power counting suggests even N4L0 
because the Pauli principle forbids three body forces without derivatives. However, such 
naive counting was already fallacious in the two body sector because of the presence of 
a low lying bound state and of a linear divergence in each of the bubble graphs making 
up Fig. |l|. We therefore investigate more carefully the UV behaviour of the amplitude 
A{k,p) at half off-shell momenta p ^ k,'-f. The integral equation simplifies then to 
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where A(s) = — | in the spin quartet channel, and 1 in the doublet. Qi is the Legendre 
Polynomial of the second kind, / the angular momentum of the partial wave investigated. 
The equation is easily solved by a Mellin transformation, ^(^ j,)(0,p) = p-i+'^oC'**)^ and 
indeed in most channels sq is real so that only one solution exists which vanishes at infinite 
momentum and hence is cut-off independent. However, the fact that the kernel of (^ is 
not compact makes one suspicious whether this is always the case. Indeed, there exists 
one and only one partial wave in which two linearly independent solutions are found: the 
doublet S wave (triton) channel, where so{l = 0, s = |) = ±1.0062 . . . i. Therefore in this 
channel, any superposition with an arbitrary phase S is also a solution, and we find that 



A=o,s=i(^-^0,p) oc 



cos[1.00621np + 5] 
p 



(6) 



Because each value of the phase provides a different boundary condition for the solution 
of the full integral equation Fig. ^, the on-shell amplitude A{k,p = k) will also depend 
crucially on 6. This means that the on-shell amplitude seems sensitive to off-shell physics. 



and what is more, to a phase which steins from arbitrarily high momenta. A numerical 
study of the full half off-shell amplitude confirms these findings, Fig. ^. This cannot be. 
Because physics must be independent of the cut-off chosen, this sensitivity of the on- 
shell amplitude on UV properties of the solution to (|^) must be remedied by adding a 
counter term. And since the power counting in the two nucleon sector is fixed, one can 
show that a necessary and sufficient condition to render cut-off independent results is to 
promote the three body force (|) to LO, H{6) ~ Q~^, and absorb all phase dependence 
into it, see Fig. ^ The variation of H{6) with the cut-off or phase is known analytically 
from d^), but its initial value is unknown. Therefore, one physical scale S must be de- 
termined experimentally. This one, new free parameter explains why potential models 
which provide an accurate description of NN scattering can vary significantly in their 
predictions of the triton binding energy B^ and three body scattering length Og in the 
triton channel, although all of them lie on a curve in the (S3, % ) plane, known as the 
Phillips line. Fig. |^. Determining H{S) by fixing the three body scattering length to its 
physical value, the triton binding energy is found in ENT(^) to be 8.0 MeV at LO. 
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Figure 6. Left: Variation of the cut-off by a factor 3 in a numerical study of the naive 
off-shell amplitude Fig. ^ in the triton channel has a dramatic effect on the scattering 
length a{p = 0). Right: Comparison of the ENT prediction for the Phillips line with 
results from various potential models. Figures from Ref. [|1^ . 
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Figure 7. The cut-off independent Faddeev equation in the triton channel. 



It must again be stressed that the three body force of strength H{6) was added not out 
of phenomenological needs. It cures the arbitrariness in the off-shell and UV behaviour 
of the two body interactions which would otherwise contaminate the on-shell amplitude. 
Just as the off-shell behaviour of the two body amplitude Fig. ^, the strength H{6) is 
arbitrary, and only the sum of two and three body graphs is physically meaningful. 



Summing the deuteron bubbles, each graph in the upper hne of Fig. ^ behaves hke p~^ 
in the UV. But the solution of the Faddeev equation goes like p~^+^o with irrational (or 
even complex) sq{1, s) [0 and is hence more than the naive sum of graphs. 

The limit cycle thus encountered in the triton channel is a new renormalisation group 
phenomenon and also explains the Efimov and Thomas effects p|JT6|l. However, the phe- 
nomenon occurs only in the triton channel. In all other partial waves, three body forces 



enter only at higher orders [T^ 



2.3. Partial Waves in Neutron— Deuteron Scattering 

In the three body sector, the equations to be solved in ENT(]f) and ENT(KSW) are 
computationally trivial and can furthermore be improved systematically by higher order 
correction which involve only (partially analytic, partially numerical) integrations, in 
contradistinction to many- dimensional integral equations arising in other approaches. A 
comparative study between the theory with explicit, perturbative pions (ENT(KSW)) 
and the one with pions integrated out was performed [Q in the spin quartet S wave for 
momenta of up to 300 MeV in the centre-of-mass frame (-Ecm ~ 70 MeV). As seen above, 
the two formulations are identical at LO. Because three body forces enter only at high 
orders, this channel is completely determined by two body properties at the first few 
orders and no new, free parameters enter. 

The calculation with/ without explicit pions to NL0/N2L0 shows convergence: For 
example, the scattering length is a(^S'3 , LO) = (5.1 ± 1.5) fm, and at NLO with (without) 
perturbative pions a(^^3, NLO, KSW)'= (6.8±0.7) fm (a(4^3, NLO,^^) = (6.7±0.7) fm). 
At N2L0, [|16| report 0(^6*3, N2L0, |f) = (6.33 ± 0.1) fm, and the experimental value is 
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5*3 , exp) = (6.35±0.02) fm. Comparing the NLO correction to the LO scattering length 

provides one with the familiar error estimate at NLO: (|)^ ~ 10%. The NLO calculations 
with and without pions lie within each other's error bar. The N2L0 calculation is inside 
the error ascertained to the NLO calculation and carries itself a theoretical uncertainty 
of about (|)^ ~ 4%. The calculation of pionic corrections in ENT(KSW) shows that they 
are - although formally NLO - indeed much weaker. The difference to ENT(^) should 
appear for momenta of the order of nij^ and higher because of non-analytic contributions 
of the pion cut, but those seem to be very moderate, see Fig. |^. This and the lack of data 
makes it difficult to assess whether the KSW power counting scheme to include pions as 
perturbative increases the range of validity over the pion-less theory. 
Finally, the real and imaginary parts of the higher partial waves Z = 1, ... ,4 in the spin 
quartet and doublet channel were found f^ in a parameter-free calculation, see Fig. ^j. 
Within the range of validity of this pion-less theory, convergence is good, and the results 
agree with potential model calculations (as available) within the theoretical uncertainty. 
That makes one optimistic about carrying out higher order calculations of problematic 
spin observables like the nucleon-deuteron vector analysing power Ay where EFT will 
differ from potential model calculations due to the inclusion of three-body forces. 

3. Outlook 

Many questions remain open: Does the power counting Nature chose include pions 
perturbatively or non-perturbatively? At the moment, some people advocate a mixture. 
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Figure 8. Real parts of the quartet S 0] and doublet D wave phase shifts in nd 
scattering. Legend left: Dashed: LO; solid (dot-dashed) line: NLO with perturbative 
pions (pions integrated out); dotted: N2L0 without pions. Realistic potential models: 
squares, crosses, triangles. Stars: pd phase shift analysis. 



To investigate processes with pions in the initial or final state might be helpful [|T8|. How 
to extend the analysis in the triton channel systematically to higher orders? Technically, 
how to regularise a Faddeev equation numerically with external currents coupled in a field 
theory? How do four body forces scale? Extend to nuclear matter! 
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